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The bulk viscosity, £ and its ratio with the shear viscosity, Q/r) have been studied in an anisotrop- 
ically expanding pure glue plasma in the presence of turbulent color fields. It has been shown that 
the anisotropy in the momentum distribution function of gluons, which has been determined from a 
linearized transport equation eventually leads to the bulk viscosity. For the isotropic (equilibrium) 
state, a recently proposed quasi-particle model of pure SU(3) lattice QCD equation of state has 
been employed where the interactions are encoded in the effective fugacity. It has been argued that 
the interactions present in the equation of state, significantly contribute to the bulk viscosity. Its 
ratio with the shear viscosity is significant even at 1.5T C . Thus, one needs to take in account the 
effects of the bulk viscosity while studying the hydrodynamic expansion of QGP in RHIC and LHC. 
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I. INTRODUCTION 

It is by now well established that Quark-gluon plasma 
(QGP) has been created in RHIC experiments, and 
is a strongly coupled fluid [l(. There have been first 
few reports of QGP in Pb-Pb collisions @2.76 Tev in 



LHC[6|, which reconfirm the formation of strongly cou- 
pled fluid. QGP at RHIC has shown a robust collective 
phenomenon, viz., the elliptic flow[2j. In the heavy-ion 
collisions at LHC, there are other interesting flows, viz., 
the dipolar, and the triangular flow which are sensitive to 
the initial collision geometry |3j. In this concern, we re- 
fer the reader to the very recent interesting studies [4J, |5( , 
where these new kind of flows at LHC have been investi- 
gated. 

The shear and bulk viscosities (jj and £) characterize 
dissipative processes in the hydrodynamic evolution of a 
fluid. The former accounts for the entropy production 
due to the transformation of the shape of hydrodynamic 
system at a constant volume. On the other hand, lat- 
ter accounts for the entropy production at the constant 
rate of change of the volume of the system (in the con- 
text of RHIC the system stands for the fireball). These 
transport parameters serve as the inputs from the hydro- 
dynamic evolution of the fluid. Their determination has 
to be done separately from a microscopic theory (either 
from a transport equation with appropriate force, colli- 
sion and source terms or from the field theoretic approach 
using Grccn-Kubo formula). It has been found that QGP 
possess a very tiny value of the shear viscosity to entropy 
density ratio, rj/s pj. On the other hand, bulk viscos- 
ity has achieved considerable attention in the context of 
QGP in RHIC after the interesting reports on its rising 
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value close to the QCD transition temperature [3,Q. In 
the recent investigations, these transport coefficients are 
found to be sensitive to the interactions [HJ, LLl| , and na- 
ture of the phase transition in QCD 12 1. 

The computation of transport coefficients in lattice 
QCD is a very non-trivial exercise, due to several un- 
certainties and inadequacy in their determination. De- 
spite, there are a few first results computed from lattice 
QCD for bulk and shear viscosities [13, LL6( which have 
observed a small value of rj/s, and a large value for £/ s 
at RHIC. While determining the behavior of the spectral 
function in [13| , a contribution coming from a ^-function 
has not been taken in to account. This issue has been 
discussed extensively in [1J|. The spectral density has 
been modified by incorporating the contributions from 
the (5-function by Meyer in |15| . However, a more refined 
lattice studies on 77 and £ are awaited in the near future 
with less dependence on the lattice artifacts and uncer- 
tainties. Subsequently, the possible impact of the large 
bulk viscosity of QGP in RHIC have been studied by sev- 
eral authors; Song and Heinz [17| have studied, in detail, 
the interplay of shear and bulk viscosities in the context 
of collective flow in heavy ion collisions. Their study re- 
vealed that one can not simply ignore the bulk viscosity 
while modeling QGP in heavy ion collisions. In this con- 
text, there are other interesting studies reported in the 
literature [18l - t24| . The role of bulk viscosity in freeze out 
phenomenon has been reported in [20Ll25| . Effects of bulk 
viscosity in hadronic phase, and in the hadron emission 
have been reported in [2g|. There has been a wealth of 
recent literature on the computations of bulk viscosity in 
the context of cosmology [27| , strange quark matter [28| , 
and neutron stars [29[ ■ 

The noteworthy point is that most of works devoted 
to study the hydrodynamic evolution of QGP, employ 
constant value of rj/s [3(| an d C/ s [2l|- This may not be 



desirable, in the light of experimental and phcnomcnolog- 
ical observation for QGP at RHIC. The work presented in 
this paper is an attempt to achieve, (i) temperature de- 
pendence of transport coefficients, in particular, £, (ii) to 
understand the large bulk viscosity of QGP. In this study, 
we shall take inputs from the com put ations of bulk vis- 
cosity in quasi-particle models [32l . |33[ , and combine the 
understanding with a transport theory determination of 
C in the presence of Chromo-Weibel instabilities [34|, biff ■ 
In this context the shear viscosity of QGP has already 
been addressed [T(J, [H 0, [35| , and we find very interest- 
ing results. As it is well emphasized by Pratt [37| that 
there may be a variety of physical phenomena which can 
lead to viscous effects in QGP. Among them, in this pa- 
per, we are particularly interested in the viscous effects 
which get contributions from the classical chromo-ficlds. 

The idea adopted here is based on the mechanism, 
earlier proposed to explain the small viscosity of a 
weakly coupled, but expanding hot QCD plasma J34ll35|. 
This mechanism is based on the particle transport the- 
ory in turbulent plasmas [38[ which arc characterized 
by strongly excited random field modes in the cer- 
tain regimes of instability, which coherently scatter the 
charged particles and thus reduce the rate of momentum 
transport. This eventually leads to the suppression of the 
transport coefficients in plasmas. This phenomenon in 
electro-magnetic (EM) plasmas has been studied in j39| . 
and generalized by Asakawa, Bass and Muller [3J] to the 
Non-Abelian plasma (QCD), and further employed for 
the realistic QGP EOS in [13, [HI- As it is emphasized 
in [401 ! tne sufficient condition for the spontaneous forma- 
tion of turbulent, partially coherent fields is the presence 
of instabilities in the gauge fields due to the presence of 
charged particles. This condition is met in both EM plas- 
mas with an anisotropic momentum distribution [4l| of 
charged particles and in QGP with an anisotropic distri- 
bution of thermal partons [42[ • Here, we shall argue that 
the similar mechanism can lead to a large bulk viscosity 
for the hot QCD plasma for the temperatures relevant at 
RHIC and heavy ion collisions at LHC. 

The paper is organized as follows. In Sec. II, we 
present the general formalism to determine the trans- 
port parameters from a transport equation with a Vlasov 
term. We have neglected the collision and source term, 
while obtaining bulk viscosity. In Sec. Ill, we discuss the 
temperature dependence of bulk viscosity and its com- 
parison with the shear viscosity. Finally, in Sec. IV, we 
present the conclusions and outlook. 



II. TRANSPORT PARAMETERS WITHIN A 
QUASI-PARTICLE MODEL 

The determination of transport coefficients requires 
modeling beyond the equilibrium properties, in terms 
of the collision terms and other transport parameters, 
and also the nature of perturbation to the equilibrium 
distribution. In particular, their determination within 



linearized transport theory needs knowledge of EOS and 
the equilibrium momentum distribution functions of par- 
ticles, which constitute the plasma. We shall first discuss 
the modeling of the EOS within a quasi-particle model. 
The EOS chosen here is the pure SU(3) gauge theory 
EOS [43| . We subsequently discuss the setting up of the 
transport equation and the determination of £. 



A. The quasi-particle model 

Lattice QCD is the best, and most powerful technique 
to extract non-perturbative information on the equation 
of state for QGP [44|, [45| ■ Recently, we have proposed a 
quasi-particle model to describe the lattice data on pure 
SU(3) gauge theory pressure (LEOS), and studied the 
bulk and transport properties of QGP [ljj , which is uti- 
lized in obtaining the temperature dependence of bulk 
viscosity here. In this description, quasi-gluon distribu- 
tion function extracted from LEOS possess the following 
form, 



z g exp(-ftp) 
1 - z g exp(-/3p) 



(1) 



It has further been argued [ll| that the model is in the 
spirit of Landau theory of Fermi liquids. The connec- 
tion with the Landau's theory is apparent from the single 
quasi-gluon energy, which gets non-trivial contributions 
from the quasi-particle excitations. The dispersion rela- 
tion (single particle energy) came out to be, 



E p =p + T 2 d T Hz g ). 



(2) 



The main feature of the description is the mapping 
of strongly interacting LEOS in to a system of non- 
interacting/weakly interacting quasi-gluons (free up to 
the temperature dependent fugacity, z g which encodes all 
the interactions, and the dispersion relation in Eq.(|2"])). 
This enables us to tackle highly non-trivial strong inter- 
action in QCD in a very simplified manner while study- 
ing the properties of QGP. Interestingly, Eq.®, which 
is obtained from the thermodynamic definition of the 
energy-density in terms of Grand-canonical QCD par- 
tition functions, ensures the thermodynamic consistency 
in hot QCD, and reproduces the lattice results on the 
trace anomaly correctly. This is also true for the re- 
cently proposed quasi-particle model which describes the 
(2 + l)-flavor lattice QCD g|. 

This quasi-particle understanding of hot QCD has been 
quite successful in describing the realistic QGP equa- 
tions of state, and in investi gat ing the bulk and trans- 
port properties of QGP [13 111 El 51 • We shall uti- 
lize Eqs.([T|) and @ to determine the bulk and shear vis- 
cosities within transport theory framework here. Note 
that there arc other quasi-particle approaches to describe 



lattice QCD EOS based on effective thermal masses 
for quasi-partons |49l - [53| , approaches based on Polyakov 
loop 154 . and quasi-particle models with gluon conden- 
sate [55|, [56|. Recently, transport coefficients for QGP 
within the effective mass models in the relaxation time 
approximation have been reported in [57|, |58j . As argued 
in J48[, our model is distinct from all these approaches, 
but equally successful in describing the thermodynamics 
of QGP. 



B. Determination of the transport coefficients 

We now consider the important physical quantities, the 
bulk viscosity, £, its ratio with entropy density, £/ s. For 
the entropy density, we again utilize the lattice results 



quoted in 11]. These quantities are very crucial to un- 
derstand the QGP in RHIC. Their determination requires 
knowledge of the collisional properties of the medium 
when it is perturbed away from equilibrium. To deter- 
mine these quantities, we adopt approach of Illil34 ,l35l|. 
The shear viscosity had been determined in [11], which 
we shall utilize to study the ratio ^/j] in the later part. 
Here, we consider £ and determine it from a transport 
equation. 

The determination of bulk viscosity has been done in 
a multi-fold way. Firstly, we need an appropriate mod- 
eling of distribution function for the equilibrium state. 
Secondly, we need to set up an appropriate transport 
equation to determine the form of the perturbation to the 
distribution function. These two steps eventually deter- 
mine the bulk viscosity. For the former step, we employ 
the quasi-particle model for LEOS discussed earlier. We 
shall leave the analysis in the case of full QCD for future 
investigations. 

The bulk viscosity has two contributions same as the 
shear viscosity in [341 ] . (i) from the Vlasov term which 
captures the long range component of the interactions, 
and (ii) the collision term, which models the short range 
component of the interaction. Here, we shall only concen- 
trate on the former case. The determination of shear and 
bulk viscosities from an appropriate collision term will be 
a matter of future investigations. Importantly, the analy- 
sis adopted here is based on weak coupling limit in QCD, 
therefore, the results are shown beyond 1.3T C assuming 
the validity of weak coupling results for QGP there. 



1. Formalism 

Let us first briefly outline the standard procedure 
of determining transport coefficients in transport the- 
ory [3J, [5!|. The bulk and shear viscosities, £ and r\ 
of QGP in terms of equilibrium parton distribution func- 
tions are obtained by comparing the microscopic defini- 
tion of the stress tensor with the macroscopic definition 
of the viscous stress tensor. The microscopic definition 
of the stress tensor is, 



Ti, 



d 3 P ,,„ - 

—^rPiPkf(P,r)- 



{2n) 3 E, 



(3) 



On the other hand, macroscopic expression for the vis- 
cous stress tensor is given by, 



Tk = P8 ik + eu t u k - 2ry(Vu). tfc - C^ifcV 



(4) 



where (Vm)^ is the traceless, symmetrized velocity 
gradient, and V • u is the divergence of the fluid velocity 
field. E p accounts for the dispersion relation. To deter- 
mine £ an 77, one writes the gluon distribution function 
as 



f(p,r) 



1 



z g 1 cxp(/3w -p- fi(p,r)) -1 



(5) 



Assuming that f\ (p, r) is a small perturbation to the 
equilibrium distribution, we expand f(p, r) and keep the 
linear order term in f\\ this leads to, 



= /o(p)fl + /i(p>1(l + /o(p))), 



where /o(p) is the isotropic distribution function, as we 
shall see that this will be same as the equilibrium ther- 
mal distribution function of the quasi-gluons, in the rest 
frame of the fluid. As discussed in [ll| , £ and 77 are deter- 
mined by taking the following form of the perturbation 
/1, 



flip, 



E p T 2 ' 



PiPj Ai(p)Vu)^ + A 2 (p)(V.u)^ 



(7) 

where the dimensionless functions Ai(p),A2(p) mea- 
sure the deviation from the equilibrium configuration. 
Ai(p), A2(p), lead to r\ and ( respectively. Note that 
Ai(p) is an isotropic function of the momentum in con- 
trast to A2(p), which is an anisotropic in momentum p. 
Since £ and rj are Lorentz scalars; they may be eval- 
uated conveniently in the local rest frame. In the lo- 
cal rest frame of the fluid /o = f eq . Considering the 
a boost invariant longitudinal flow, V • u — — and, 
(Vu)r/ = 4-diag(— 1, — 1, 2), in the local rest frame, we 
find that fi(p) takes the form, 



hip) 



E p T 2 t 



3 



^ v 2 

l' • 



E p T 2 t 



(8) 



where t is the proper time(r = \/t 2 — z 2 ). The shear 
and bulk viscosities are obtained in terms of entirely un- 
known function Ai(p) and A-zip) as, 



c 



w 9 

3T 2 



15T 2 



d 3 p p 2 
8^Ej 



d 3 p p 4 
8^El 



Al(p)/e,(l + /e,), 



(9) 



(p 2 -3 C ^)A 2 (p)/ eg (l + / eg ). (10) 



In these expressions, v g = 2(N 2 — 1) is the degrees 
of freedom. Notice that while obtaining the expression 
for the bulk viscosity, we have exploited the Landau- 
Lifshitz condition for the stress energy tensor. The factor 
— (3c 2 -E 2 ) in the rhs of Eq. (fTU|) is coming only because of 
that. For details, we refer the reader to [33| . 

The determinations of Ai(p), and r\ have already been 
done in [10|, [Dj • We shall utilize these results to fix the 
temperature dependence of £ in the later part of the anal- 
ysis. Now, we shall focus on the determination of the 
unknown function A.2(p) and £. 



2. Determination of A 2 (p) 

For simplicity, we consider the purely chromo-magnctic 
plasma for our analysis. The modeling of transport equa- 
tion for full chromo-electromagnetic plasma is straight 
forwardpjj, [ll], |34[ and differs by simple factors. Here, 
we only quote the mathematical form of the drift term 
and the Vlasov term (For details see [Tfl [34| ) . 

The drift term in the transport equation for the full 
chromo-electromagnetic plasma for LEOS is obtained as, 



L 2 = -(pxd p ) 2 + (pxd ff ) 



■(L 



P n2 



(L p z 



(13) 



Since L 2 contains angular momentum operator i p , 
therefore it gives non- vanishing contribution while op- 
erating on an anisotropic function of p. It will always 
lead to the vanishing contribution while operating on an 
isotropic function of p. Therefore, V^/ eg = 0. Now, we 
write the transport equation containing only those terms 
which contribute to bulk viscosity £ as, 



(7^-c 2 s )^(V-u)f eq (l + f eq ) 



9 2 C 2 



3(N 2 - l)E 2 



<B 2 >r m L 2 /i(P^)/e,(l + /e,). 

(14) 



Substituting for /1 in term of the unknown function 
A2(p) and rearranging above equation, we obtain a dif- 
ferential equation for A 2 (p) as, 



2(N 2 - l)TEl „ 2 „ „ 

Now, using the fact that L 2 only operates on the 
anisotropic function of p, we can write, 



(vd)f eq (p)=f eq (l + f eq ) 



PiPj 

E P T 



(Vu)i 



•d 



'rn-^P 



3T 2 d£/dT 



+ ( 



3^ 



af(V-u) 



A 2 (p) = 



2(N 2 -l)TE 2 p 



N c g 2 <B 2 > r m p 2 v 3 



(i--ciE;)xg(p), (16) 



(11) 



where g(p) can be determined from the following con- 
dition, 



where c 2 is the speed of sound. The other notations 
are kept same as in [111 ] . Note that < E 2 > stands 
for the chromo-electric field, r e ; relaxation time associ- 
ated with the instability J34J. In Eq. pTj) , first term con- 
tributes to the shear viscosity, second term contributes to 
the thermal conductivity, and the third term contributes 
to the bulk viscosity. Since, we are considering the purely 
chromo-magnetic plasma, so the second term will not be 
present. 

On the other hand the force term (Vlasov term) which 
we denote as V^, is obtained as[ll|,[3J] follows, 



9 2 C 2 



2(N 2 l)E 2 



<B 2 > 



(12) 



where C 2 = N c , < B 2 > denotes chromo-magnctic 
field, T m is the times scale associated with instability in 
the field, and the operator L 2 is 



l 2 g(p) = 1, 



which leads to, 



9(P) = 7; H — r^ = ln (— )• 
2 P a Po 



(17) 



(18) 



Since, at high temperature average value of the energy 
is 3 T. Employing equipartition theorem for relativistic 
massless gas, we obtain p\ = 6 T 2 . Substituting Ea. (fT8|) 
in Eq. p^l) . we obtain, 



2(N 2 - 1)TE 2 p 2 , „ PT 

A ^ = W^nTZ-^^ ~ c El) ln(^Xl9) 
N c g z < B l > T m p z 3 Po 

The determination of bulk viscosity is incomplete un- 
less we know not only the temperature dependence of 



the speed of sound square, c 2 , and the the collective con- 
tributions of quasi-particle to the single particle energy, 
T 2 c*t hi(z s ) but also the quantity g 2 < B 2 > r m . 

We determine first two quantities using the quasi- 
particle model. As from Ref.[ll[, the trace anomaly in 
terms of effective quasi-particle number density and ef- 
fective gluon fugacity reads, 



(e - 3P) _ M 



T A 



= j§{Td T ln(z g )}. 



(20) 



The thermodynamic quantities can be obtained using 
the well known thermodynamic relations. In particular, 
the energy density and the entropy density was shown to 
be in almost perfect agreement with the lattice data nj . 
We determine, c 2 by employing a method reported in 60J . 
The temperature dependence is shown in Fig. 1. 

To relate the denominator of Eq. (|19l) to the gluon 
quenching parameter, q we go the light cone frame. In 
this frame, Eq. (|19p can be rewritten as, 



A 2 (0) 



4(iV 2 - l)TE 2 p 2 

N c g 2 <E 2 + B 2 >T m p 2[ 3 s p) W 

(21) 
The gluon quenching parameter, q is related with the 
denominator of rhs of the above equation as [35j , 



q = 



2g 2 N c 

3(N 2 1) 



< E 2 + B 2 > t„ 



(22) 



Now, employing Ea. ([2"TT) in Eq.(8), we obtain the C as, 



{N 2 C - 1) 
ZT-K 2 q J 



OO /-OC 



Prdp T dp z {- 



cl E 2 ) 2 x 



ln(—)x/e,(l + /«*). 

Po 



(23) 



On the other hand, if we employ the results of [Tl| for 
Ai(p) in Eq.(10) for 77, we obtain, 



T 6 64(7V 2 - 1) 
V = — 0I2 p olyLog[6, z g ] , (24) 



Q 



3tt 2 



where N c = 3 and PolyLog[6, z g ] = J2T=i ~W~ 

Now scaling, all the quantities in the integrand in 
Eg. (|23)l by T, and rewriting Ea. (j24"|) in the form given 
below, we obtain, 



C = Kh{T/T c )- r, = Tl h (T/ Tc ), (25) 

q q 

where Ji(T/T c ), is evaluated by integrating the rhs of 
Eqs.flU]) numerically, and h{T/T c ) = ^PolyLog[6,z g }. 
The T/T c scaling of these quantities is coming from the 
temperature dependence of the effective gluon fugacity, 



§ 




FIG. 1. (Color online) The ratio of bulk viscosity, ( to the 
shear viscosity, rj as a function of temperature. The leading 
order (LO) result of C/v has been obtained from the data 
taken from Refs.[62. [63|. and shown as dashed line. For the 
sake of comparison, we have multiplied the leading order (,/r) 
by a factor of hundred. 



z g . Here, T c is taken to be 0.27 GeV [6l|. Clearly, the 
quantity which can be determined unambiguously in our 
approach is the ratio C/v = h(T lTc)lh(T jT c ). 

In the recent past, Chen et. al [62j, [63j have computed 
the leading order shear and bulk viscosities for purely glu- 
onic plasma. This is nothing but the collisional contribu- 
tion to these transport parameters for a gluonic plasma. 
It is to be instructive to compare the results on Q/rj ob- 
tained in the present work with those reported in [621 ] . 
This has been shown in Fig. 1, where both the results on 
e/r? are plotted as a function of T/T c . Note that while 
obtaining the temperature dependence of the ratio C/^j 
we have employed the two-loop expression for the run- 
ning coupling constant at finite temperature quoted in 
[62j. Quantitatively the ratio (/r/ is much smaller than 
what we have obtained from the diffusive Vlasov term. If 
we compare the two curves on the ratio f^jr] shown in Fig. 
1, wc find that in contrast to our prediction on C/Vi t ne 
leading order result suggests the near conformal picture 
of hot QCD even at lower temperatures. 

Next, we discuss the interplay of the two contributions 
to the bulk viscosity, viz., the anomalous, and the leading 
order (collisional). As it is emphasized in (34J], these two 
contributions for rj are inverse additive. Their inverse 
additivity has been argued from the additivity of vari- 
ous rates in the hot QCD medium. In the case of weak 
coupling, the former is predominant. It seems that a sim- 
ilar additivity of the inverse of two contributions to Q,viz. 
(denoted as C, a and £ c respectively) may perhaps be valid. 
This could be understood as follows: since ( a is inversely 
proportional to the q (transport rate), on the other hand 
collisional Cc will be inversely proportional to the colli- 
sion rate. Following the argument previously mentioned, 
one may write, Cy 1 = CT 1 + CT 1 ; where £r denotes the 
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FIG. 2. (Color online) Comparison of the ratio (,/ri with 
the perturbative QCD, and strongly coupled theories. The 
quantities R per t, and R s tr are defined in terms of the ratios 



and 



respectively. Here, pert stands for per- 



turbative QCD, and str stands for the strongly coupled near 
conformal theories. 



total bulk viscosity. This inverse additivity of Q ai and 
C, c at weak coupling, suggests that the collisional bulk 
viscosity (leading order) will dominate over the anoma- 
lous one, since the former is quantitatively much smaller 
than the latter. However, one has be very cautious while 
comparing these two contributions for the temperature 
ranges relevant for QGP at RHIC. This is because of the 
strongly coupled fluid like picture of QGP. At this mo- 
ment, we do not know whether the inverse additivity of 
£ will be followed at the temperatures which are closer 
to T c or not. This is a very crucial issue, and will require 
much deeper investigations, which is beyond the scope of 
the present work. Henceforth, we denote the anomalous 
bulk viscosity as ( dropping the subscript, a. 

We now proceed to discuss the temperature depen- 
dence of C/77 and £/s. 



C. Temperature dependence of C/77 and £/s 

In our analysis, determination of the temperature de- 
pendence of £ and r\ is incomplete, without the knowl- 
edge of the temperature dependence of q in QGP. This is- 
sue was addressed by fixing the temperature dependence 
of q by calculating the soft part of the energy density 
and the relaxation time associated with the instability of 
chromo-fields PJJ. To do that we take inputs from the 
phenomcnological values of q, which is known at a partic- 
ular temperature [64|]- Here, we have utilized the same 
transport equation and quasi-particle model developed 
for pure SU(3) lattice QCD EOS, as in [u|. Therefore, 
we employ the temperature dependence of rj/s to obtain 
the temperature dependence of C,/ s. This is quite easier 



FIG. 3. (Color online) C,/s and r//s as a function of tempera- 
ture. The dashed (green) line denotes £/s and solid(red) line 
denotes the r//s. 



to do, since the ratio, (/rj can easily be obtained from 
Eq. (fJSl). 

The temperature dependence of Q/rj is shown in Fig. 
1, C/v relative to perturbative QCD prediction [65|. and 
strongly coupled near conformal gauge theories |66f is 
shown in Fig. 2. On the other hand, £/s and rj/s are 
shown together in Fig. 3. Let us discuss their behavior 
one by one. From Fig. 1, it is clear that C,/r] is equally 
significant while studying the hydrodynamic evolution of 
hot QCD matter until we reach T = 2T C . As we go 
to the higher temperatures the ratio further decreases 



and eventually vanishes when c s 
relation E„ 



and the dispersion 



— p. Quantitatively, C/v ~ 2.3 at 1.3T C ; 1.0 
at 1.5T C ; 0.2 at 2.0T C . Therefore, for T > 2.5T C , one can 
ignore £ over 77. In other words, the hot QCD becomes 
almost conformal there. 

The ratio, Q/rj decreases as we increase the temper- 
ature. The decrease is quite steeper until we reach 
T = 2.0T C . For higher values of T it is much slower. It is 
hard to make clear cut statement in regard to the behav- 
ior of £//? with temperature, since by looking at Eqs. (23) 
and (24), it is clear that the behavior of C,jr\ as a func- 
tion of temperature is mainly governed by the tempera- 
ture dependence of trace anomaly (through quasi-gluon 
dispersion relation), speed of sound, c 2 and temperature 
dependence of z g and gluon quenching parameter, q. 

To compare the perturbative QCD prediction of the 
ratio e/77, we consider R per t = , 2 _ii2 1 where {c 2 s — |) 
can be thought of as the measure of conformal symmetry, 
which we call conformal measure. For scalar field theo- 
ries, Q/tj — 15(Cg — -|) 2 [63, and this has been found to be 
true for a photon gas coupled with hot matter by Wein- 
berg [68|] . The pre-factor 15 is not fixed for perturbative 
QCD but the scaling (/v ~ (c, - |) 2 is valid[65J]). Note 
that in certain strongly coupled near conformal theories 
with gravity dual the ratio ^/i] shows linear dependence 



on the conformal measure [66[ . To compare with the lat- 
ter, we consider the ratio R s t r = , 2 _i, ■ We have shown 

V\ c s 3 ) 

the behavior of R per t and R s t r as a function of temper- 
ature in Fig. 2. Clearly, none of these two scaling are 
respected by the ratio C/r] in Fig. 2 even at 2.5 T c . It 
is safer to say that (/rj for LEOS which is obtained from 
transport equation with Vlasov-Dupree term [n], |34[ nei- 
ther shows linear nor the quadratic dependence with the 
conformal measure, (c~ — |). However, one can realize 
the quadratic scaling of £/rj with the conformal measure 
in a certain limiting case. It is easy to say from Eqs.(24) 
and (25) that for E p = p (p << T 2 d T {ln(z g ))), if the 
thermal distribution of quasi-gluons shows near ideal be- 
havior, and with constant value of q/T , the quadratic 
scaling can be achieved. Moreover, this may perhaps be 
realized at higher temperatures which are not relevant 
for QGP in RHIC and LHC. If we compare qualitatively 
our prediction of R per t , and R s t r with the leading order 
result of [62J (see Fig. 4 of this Ref.) , we find oppo- 
site trend of these quantities at very high temperature. 
The former decreases, although slowly, in contrast to the 
latter, as a function of temperature. This could perhaps 
be because of their origin from the distinct physical pro- 
cesses in hot QCD medium. The slow decreases of the 
former at higher temperatures, could be understood as 
the effect of thermal distribution function of quasi-gluons 
(through z g , since z g increases very slowly as a function of 
temperature, and will asymptotically approach to unity) . 
Finally, in Fig. 3, we have shown the temperature 
dependence of C,/ s and rj/s. The C,/ s decreases as with 
increasing temperature for T > 1.5, in contrast to rj/s. 
As mentioned earlier, £/s and rj/s becomes equal around 
1.5T C (below which £/s is higher, and lower for higher 
temperatures.). Again the behavior is predominantly 
controlled by the behavior of c 2 s , and the trace anomaly 
through the modified dispersion relation with tempera- 
ture. 



III. CONCLUSIONS AND FUTURE 
PROSPECTS 

In conclusion, we have estimated the temperature de- 
pendence of bulk viscosity to entropy density ratio (C/s), 
and bulk viscosity relative to shear viscosity, (/r/ within 
a quasi-particle model for pure glue QCD at high tem- 
perature by employing transport theory. We have deter- 
mined C/Vi exactly and unambiguously. In our analysis, 
these quantities get contributions from the instabilities in 
the chromo-electromagnetic fields due to the anisotropic 
thermal distribution of the partons in QGP. The mecha- 



nism has succeeded in explaining the small rj/s and large 
value of the ratio C^/rj. In fact, (/tj is around 2.3 at 1.3 T c , 
of the order of unity at 1.5 T c , and 0.2 at 2 T c . This tells 
us that the breaking of conformal symmetry in hot QCD 
plays crucial role even at 2 T c . In consequence, shear 
and bulk viscosities are equally important while studying 
the hydrodynamic evolution of QGP at RHIC and LHC. 
One cannot simply ignore bulk viscosity even at 2.0T C 
while modeling the heavy ion collisions. Moreover, rj/s 
increases as a function of temperature, in contrast to £/s 
beyond 1.5T C . As expected £/s and C,jr\ are vanishingly 
small beyond 2.5T C . This may be due to the fact that 
conformal measure is very small there, and the speed of 
sound is closer to 1/3. We have compared our predictions 
on C/77 to the leading order result on the same quantity 
obtained by [621 ] . Interestingly, in the perturbative region 
(temperatures beyond 1.5T C ), our study also agree with 
the near conformal picture of hot QCD similar to leading 
order results of Chen et. al [62j • On the other hand the 
predictions are in contrast at lower temperatures. How- 
ever, this may not be thought of as the complete story, 
an adequate analysis on the interplay of our predictions 
on £, and leading order prediction is very much desired, 
and will be a matter of future investigations. 

We have addressed the temperature dependence of the 
bulk and shear viscosities of pure glue sector of hot QCD 
only. An extension to full QCD including collision term, 
employing the understanding of [46[, will be a matter of 
future investigations. We strongly believe that a similar 
analysis will also be valid in the case of full QCD. The 
most interesting study would be to include the temper- 
ature dependence of rj/s and (,/ s in the existing hydro 
codes to model QGP, and see how various obscrvables 
get modifications. Moreover, future directions may in- 
clude exploration on the effects of r\ and C, on the quarko- 
nia suppression in heavy ion collisions along the lines 
of [63, lZQ| . Finally, it would be of interest to include the 
baryon chemical potential utilizing the very recent lattice 
studies [2ill72], and determine the transport coefficients. 
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